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ABSTRACT 


The problem of Lambert is reduced to the solution of a single equation 
which is valid for elliptic, parabolic, and hyperbolic motion. The equation has 
no singularity for 180 ° transfers, holds for any number of revolutions, and ex- 
presses a normalized time of flight as a function of one parameter and a uni- 
versal two- body variable. 
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A UNIVERSAL SOLUTION OF LAMBERT’S PROBLEM 


I. INTRODUCTION 

Consider a particle moving in a gravitational, inverse- square, central force 
field and let f and r be respectively the true anomaly and the distance of the 
particle from the center of attraction at time t . If subscripts 1 and 2 indicate 
values at times tj and t 2 , Lambert’s problem assumes that t lr r 2 , and f 2 - fj 
are known and seeks to find some parameter of the motion such as the semi- 
major axis a which will enable us to compute directly the velocity components 
at times t x and t 2 , provided the plane of motion is known. 

The problem of Lambert is reduced to the solution of two equations in two 
unknowns by Battin (1964) and Pitkin (1968). These equations are universally 
valid for the three types of conic motion. Battin (1964) reports the work of 
Deyst, wherein the problem is reduced to one universally valid equation but with- 
out provision for more than one revolution during the flight time t 2 - t x . Pines 
(1965) has independently derived the equation of Deyst, with a modification which 
allows for any number of revolutions. The equation of Deyst and Pines expresses 
the time of flight as a function of two parameters and a universal variable. 

In this paper Lambert's problem is reduced to a single universal equation 
which is simpler than the equation of Deyst and Pines. It expresses a normalized 
time of flight as a function of only one parameter and a universal variable. 

We will first present separate equations for the elliptic and hyperbolic cases. 
We will then combine these equations into a single universal equation. 

II. ELLIPTIC AND HYPERBOLIC CASES 

For elliptic motion we have 


r - a(l - e cos E), 

(1) 

(/j/a 3 ) 1/2 (t - t p ) = E - e sin E, 

(2) 

r 1/2 sin-yf = [a(l + e )] 1/2 sin^E , 

(3) 

r^oos-jf = [a(l “ e)] 1/2 cos "j E , 

(4) 
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where ^ is a known constant, a is the semi-major axis, E is the eccentric 
anomaly, t is a time of periapsis passage, and e is the eccentricity of the 
orbit. We define 


a 2 ( E 2 " E l) ’ 

oo - e cos ~2 (Ej + E 2 ) . 


Place subscripts 1 and 2 on r and E in (1) and add the two resulting 
equations to obtain 


( f i + f 2)/ 2a = 1 " u 003 a ■ 


( 5 ) 


Place subscripts 1 and 2 on t and E in (2) and subtract the resulting 
equations to obtain 


•j (/Va 3 ) 1/2 (t 2 - t x ) = a - go sin a . (6) 

Place subscripts 1 and 2 on r , f , and E in (3) and (4) and combine the 
resulting equations to obtain 

(i-jT^ 172 cosy (f 2 - f,) = a(oos a - a>) . (7) 

If we define 

T = (M I/2 (t 2 - »,) / ( r, + r 2 ) 3/2 , 

q = [2 ( r i r j) ,/2 /( r i + r 2 )] “ s T - f i) 

and eliminate a and 00 from (5), (6), and (7), we get 
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T - (~l) m (1 ~ q cos a) 1/2 [q + (a ~ sin a cos a) (1 - q oos a)/sin 3 a] , (8) 

where m is the number of complete revolutions which occur during the time of 
flight. We note that 

-1 < q < 1 . 

An equation similar to (8) has been derived by Godal (1961) from formulas 
due to Gauss, but with the flight time expressed as a function of two parameters 
and the eccentric anomaly difference. 

Eliminating <u from (5) and (7) gives 

( r i + r 2) (1 ” d a>s a ) ~ 2a sin 2 a (9) 

for the calculation of a after a has been found from (8) . 

If q is near 1, it may be desirable to replace 1 - q cos a by 
2 sin 2 \ a + (1 - q) cos a in (8) and (9) for better numerical accuracy. 

Kepler's equation for elliptic motion can be written in the form 

~2 (^/a 3 ) 1/2 (t 2 - tj) = a + jVj fjy^/ja) 172 ^ sin 2 a - (l - sin a cos a, (10) 

where fj = dr/dt at time tj . Having found a and a, we can obtain r t from (10). 
In a similar manner we obtain for the case of hyperbolic motion 

T - (1 ~ q cosh / 3) 1/2 [q “ (J3 - sinh /3 oosh /?)( 1 - q cosh /3)/sinh 3 /3] , (H) 

where q and T have the same definitions as before and 

e = -jK -h.) , 

H being the usual hyperbolic parameter. The analogue of (9) is 
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( r i + r 2 ) (1 " Q cosh 1 3 ) - -2a sinh 2 fi. (12) 

We note that we have adhered to the sign convention for the semi-major axis a, 
i.e., a > 0 for an ellipse and a < 0 for a hyperbola. 

Kepler's equation for hyperbolic motion can be written in the form 


\ (/V-a 3 ) 1/2 (t 2 - tj) = + [( r i Tj j(-fu.a) 1 /2 j sinh 2 /? + (l “ rja) sinh/?cosh/3,(13) 

from which fj can be obtained. 

The following formulas for the eccentricity, semilatus rectum, and the 


component of velocity perpendicular to the radius vector are valid for all 

forms of two-body motion: 

e 2 

= (1 * ya ) 2 + (Vj) 2 /W. 

(14) 

P 

= 2r x - r^/a - (r^) 2 //., 

(15) 

u i 

= (fjp) 1/2 /r 1 . 

(16) 


IH. THE UNIVERSAL FORM 

We begin by defining two series 


00 

S(y) = J~ 1 (-l)V/(2i + 3)!, 

CD 

= V(-l)V/(2i + 2)! . 

(17) 

i=0 

i*0 


convergent for all values of y . With these series we can write 


sin a = a(l “ 7 S), 

00 s a = 1 - 7 C if y ~ a 2 , 

(18) 

sinh fi - /3( 1 — -y S) , 

x>sh fi - 1 ~ yC if y - - fi 2 ■ 

(19) 
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Substituting these expressions into (8) and (11) gives 


T = (“1 ) m y 1/2 [q + y (Cx + S)/x 3 ] (20) 

x = 1 - yS, y - 1 - q (1 - yC) = >C + (1 - q)( 1 - yC ) • 


Equation (20) replaces (8) and (11) if we impose the convention that y > 0 for 
elliptic motion and y < 0 for hyperbolic motion. Equation (20) becomes truly 
universal if we agree that y - 0 for parabolic motion. Of course m = 0 for 
hyperbolic and parabolic motions, y is the universal two-body variable used in 
the universal solution of the two-body initial-value problem as described by 
Battin (1964). 

Substituting (18) and (19) into (9) and (12) gives 

j = 2yx ! z , 


where 


Z = y/(r, + r 2 ) . 

Substituting (18) and (19) into (10) and (13) gives 
r j r j = yi{t 2 - tj) z - (-l) m (fi/ 2z) 1/2 [(Cx + S)/x 3 + 2t x (1 - ^Ot)] . (22) 

The following identities are useful for computing the C and S functions for 
large values of y: 


2C(4y) = [1 “ tS(7)] 2 , 

4S(4 y) = S(y) + C(y) [1 - yS(y)1 . 

Since the C and S functions are evaluated so frequently, it is worthwhile to 
develop approximations that minimize the number of arithmetical operations 
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required. Herron, et al. (1968) have developed several such approximations, 
based on theories of Chebyshev and Knuth, with bounds for the errors incurred 
when using them. 
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